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Abstract 

The energy loss of fast partons traversing the strongly interacting matter produced in 
high-energy nuclear collisions is one of the most interesting observables to probe the 
nature of the produced medium. The multipurpose Compact Muon Solenoid (CMS) 
detector is well designed to measure these hard scattering processes with its high res- 
olution calorimeters and high precision silicon tracker. Analyzing data from pp and 
PbPb collisions at a center-of-mass energy of 2.76 TeV parton energy loss is observed 
as a significant imbalance of dijet transverse momentum. To gain further understand- 
ing of the parton energy loss mechanism the redistribution of the quenched jet energy 
w^as studied using the transverse momentum balance of charged tracks projected onto 
the direction of the leading jet. In contrast to pp collisions, a large fraction the momen- 
tum balance for asymmetric jets is found to be carried by low momentum particles at 
large angular distance to the jet axis. Further, the fragmentation functions for leading 
and subleading jets were reconstructed and were found to be unmodified compared 
to measurements in pp collisions. The results yield a detailed picture of parton prop- 
agation in the hot QCD medium. 
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Figure 1: A<p-[2 distributions for leading jets of pT,i > 120 GeV/c with subleading jets of pj^2 > 
50 GeV/c for 7 TeV pp collisions (a) and 2.76 TeV PbPb collisions in several centrality bins: (b) 
50-100%, (c) 30-50%, (d) 20-30%, (e) 10-20% and (f) 0-10%. Data are shown as black points, 
while the histograms show^ (a) PYTHIA events and (b)-(f) PYTHIA events embedded into PbPb 
data. The error bars show the statistical uncertainties. 

1 Introduction 

Heavy ion collisions at the Large Hadron Collider (LHC) are expected to produce matter at en- 
ergy densities exceeding any previously explored in experiments conducted at particle acceler- 
ators. One of the key experimental signatures suggested to study the quark-gluon plasma(QGP), 
which is expected to be formed at these high energy densities, is the energy loss of high- 
transverse-momentum partons passing through the medium fT|. This parton energy loss is 
often referred to as "jet quenching". The energy lost by a parton provides fundamental infor- 
mation on the thermodynamical and transport properties of the traversed medium. 



2 Data set and experimental method 

The dijet analysis presented in this paper was performed using the data collected from pp and 
PbPb collisions at a nucleon-nucleon center-of-mass energy of ^/s~^= 2.76 TeV at the Com- 
pact Muon Solenoid (CMS) detector. The CMS detector has a solid angle acceptance of nearly 
An and is designed to measure jets and energy flow^, an ideal feature for studying heavy ion 
collisions. The CMS detector is described in detail elsewhere [2J. This analysis is based on a 
total integrated PbPb (pp) luminosity of 6.7 jib^^ (225 nb^^). The jets used for analysis were 
reconstructed with two different algorithms. The results related to dijet correlations uses jets 
reconstructed from calorimeter energies using an iterative cone algorithm with a cone radius 
of 0.5. The fragmentation function results use jets reconstructed using the CMS particle flow al- 
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Figure 2: Dijet asymmetry ratio, Aj, for leading jets of px,! > 120 GeV/c, subleading jets of 
Pt,2 >50 GeV/c and A(/)i2 > 2n/3 for 7 TeV pp collisions (a) and 2.76 TeV PbPb collisions in 
several centrality bins: (b) 50-100%, (c) 30-50%, (d) 20-30%, (e) 10-20% and (f) 0-10%. Data are 
shown as black points, while the histograms show (a) PYTHIA events and (b)-(f) PYTHIA events 
embedded into PbPb data. The error bars show^ the statistical uncertainities. 



gorithm, which combines information from all detector systems ||3l and the anti-k^ algorithm, 
as encoded in the Fastjet framework, with a resolution parameter of R = 0.3 [H- For both 
algorithms the underlying event is subtracted using the method described in [51 and the re- 
constructed jet energies are corrected for the detector response based on pythia |I6| simulations. 
Charged particles are reconstructed in the CMS Silicon tracking system ||2|. 



3 Dijet Correlations 

For this analysis dijets are selected with a leading jet of px^ > 120 GeV/c, a subleading jet 
of px,2 >50 GeV/c both within \f]\ < 2 and with an opening angle of A(pi2 > 2ti/3, unless 
otherwise noted. More details on the analysis presented in this section can be found in JT). 

3.1 Dijet azimuthal correlations 

One possible medium effect on the dijet properties is a change of the back-to-back alignment 
of the tw^o partons. This can be studied using the event-normalized differential dijet distribu- 
tion, {1 / N){dN / dA(p-[2), versus Acpu- Figure [ijshows distributions of Acpu between leading and 
subleading jets which pass the respective px selections. In Figure [l] (a), the dijet Acpu distribu- 
tions are plotted for 7 TeV pp data in comparison to the corresponding pythia simulations 
using the anti-A;x algorithm for jets based on calorimeter information. PYTHIA provides a good 
description of the experimental data, w^ith slightly larger tails seen in the PYTHIA simulations. 



3.2 Dijet momentum balance 
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Figure 3: Left Panel: Average missing transverse momentum, {^j), for tracks with px > 

0.5 GeV/c, projected onto the leading jet axis (solid circles). The {/j.) values are shown as a 

function of dijet asymmetry Aj for 0-30% central events. Middle Panel: The (p'j) values as a 

function of Aj inside (AR < 0.8) one of the leading or subleading jet cones. Right panel: (//|) 
outside (AR > 0.8) the leading and subleading jet cones. For the solid circles, vertical bars 
and brackets represent the statistical and systematic uncertainties, respectively. Colored bands 

show the contribution to (^|) for five ranges of track pj. For the individual px ranges, the 
statistical uncertainties are shown as vertical bars. 



Figures[T](b)-(f) show the dijet Acpu distributions for PbPb data in five centrality bins, compared 
to PYTHIA+DATA simulations. For all centrality bins a good agreement with the pythia+data 
reference is found for the bulk of the data, especially for A(pi2 > 2. 

3.2 Dijet momentum balance 

To characterize the dijet momentum balance (or imbalance) quantitatively, we use the asym- 
metry ratio, 

Pj^^Pj^^ (1) 

Pt,i + Pt,2 

where the subscript 1 alw^ays refers to the leading jet, so that Aj is positive by construction. In 
Fig. 0(a), the Aj dijet asymmetry observable calculated by pythia is compared to pp data at ^/s 
= 7 TeV. Data and event generator are found to be in excellent agreement. The centrality depen- 
dence of Aj for PbPb collisions can be seen in Figs. |2] (b)-(f), in comparison to PYTHIA+DATA 
simulations. The dijet momentum balance exhibits a dramatic change in shape for the most 
central collisions. In contrast, the PYTHIA simulations only exhibit a modest broadening, even 
when embedded in the highest multiplicity PbPb events. This observation is consistent with a 
degradation of the parton energy, or jet quenching, in the medium produced in central PbPb 
collisions. 

3.3 Overall momentum balance of dijet events 

Information about the overall momentum balance in the dijet events can be obtained using the 
projection of missing px of reconstructed charged tracks onto the leading jet axis. For each 
event, this projection was calculated as 



4 Fragmentation Functions 



A = J2-PJ cos (^i - ^Leading Jet ) , (2) 

i 

where the sum is over all tracks with pj > 0.5 GeV/c and \tj\ < 2.4. The results were then 

averaged over events to obtain (^j). The leading and subleading jets were required to have 

\r]\ < 1.6. The left panel of Fig. b^shows {^j) as a function of Aj for the 0-30% most central 
PbPb collisions. Using tracks with \f]\ < 2.4 and px > 0.5 GeV/c, one sees that, even for events 
with large observed dijet asymmetry, the momentum balance of the events, show^n as solid 

circles, is recovered within uncertainties. The figure also shows the contributions to {]/j.) for 
five transverse momentum ranges from 0.5-1 GeV/c to px > 8 GeV/c. The vertical bars for each 

range denote statistical uncertainties. A large negative contribution to (p'j) (i.e., in the direction 
of the leading jet) by the pj > 8 GeV/c range is balanced by the combined contributions of low^- 
momentum tracks from the 0.5-2 GeV/c regions. 

Further insight into the radial dependence of the momentum balance can be gained by studying 

(p'j) separately for tracks inside cones of size AR = 0.8 around the leading and subleading jet 
axes, and for tracks outside of these cones. The results of this study are shown in Fig. |3] for 
the in-cone balance (middle panel) and out-of-cone balance (right panel). As the underlying 
PbPb event is not (/^-symmetric on an event-by-event basis, the back-to-back requirement w^as 
tightened to Acpu > 5n/6 for this study. 

One observes that the in-cone imbalance of (p'j) ~ —20 GeV/c is found for the Aj > 0.33 

selection, which is balanced by a corresponding out-of-cone imbalance of (^|) ~ 20 GeV/c. 
In the PbPb data the out-of-cone contribution is carried almost entirely by tracks w^ith 0.5 <px 
< 4 GeV/c . 

4 Fragmentation Functions 

Fragmentation functions are measured by correlating tracks within the jet cones of recon- 
structed dijets with the jet axis for each jet. For this analysis dijets are selected with a leading jet 
of px,i > 100 GeV/c, a subleading jet of px,2 >40 GeV/c within \rj\ < 2 and an opening angle 
of A(/)i2 > In/ 3. The reconstructed jets are corrected to the generator final state particle level 
and the reconstructed charged particles are corrected for tracking efficiency. The fragmentation 
function is a function of ^, defined as: 

f=-lnz = -ln^, (3) 

Pt 

where px ^'^^ is the transverse momentum of the reconstructed jet, and px ^''"'^'^ is the transverse 
momentum of tracks w^ithin the jet cone (AR = y^Acp'^ + A?/^ = 0.3 around the jet axis). Tracks 
are selected with px ^''"'^^ > 4 GeV/c, which restricts this measurement to the low ^ region of the 
fragmentation function. For the PbPb analysis, this track selection minimizes the contribution 
of the underlying event to the jet fragmentation function. Figure |4] shows the reconstructed 
leading and subleading jet fragmentation functions (top row) in a cone of AR = 0.3 around 
the respective jet axis for pp collisions (left panel) peripheral PbPb (middle panel) and central 
PbPb (right panel). The corresponding jet px distributions are shown in the bottom row. 

For a direct comparison between pp and PbPb, the pp data has to include the resolution dete- 
rioration due to the underlying event fluctuations seen in PbPb. For this purpose the pp jet px 
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Figure 4: Fragmentation functions reconstructed in pp, peripheral PbPb, and central PbPb data, 
with the leading (open circles) and subleading (solid points) jets shown in the top row. The 
middle row is the ratio of each PbPb fragmentation function to its pp reference. The pp refer- 
ence is constructed by convoluting the pp fragmentation with jet pj resolution effects of PbPb 
collisions and reweighting with jet pj. Error bars are statistical, the hollow boxes represent the 
systematic uncertainty for the leading jet, and yellow boxes are the systematic uncertainty for 
the subleading jet. The bottom row shows the jet px distributions of the samples (pp, periph- 
eral PbPb, and central PbPb). It is important to note that jet pj distributions of leading and 
subleading jets are different. 

has been artificially smeared by the fluctuations observed in PbPb collisions taking into account 
the mean and RMS of the fluctuations as well as the correlation of the background level seen in 
back-to-back jet cones and the centrality dependence of the fluctuation parameters. The ratio 
between the PbPb and fluctuation convoluted pp reference fragmentation functions are show^n 
in the middle rows for peripheral and central PbPb collisions. In both cases, the fragmentation 
functions of the leading and subleading jet are in agreement with pp collisions. 
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Figure 5: Fragmentation functions in Aj bins, reconstructed in central PbPb and pp reference 
for the leading (open circles) and subleading (solid points) jets. The middle row shows the 
ratio of each fragmentation function to its convoluted pp reference. The systematic uncertainty 
is represented by an hollow box (leading jet) or yellow box (subleading jet). Error bars shown 
are statistical. The bottom row shows the jet px distributions. The reference distributions are 
reweighted to match the jet pj distributions in PbPb. 

4.1 Fragmentation Functions in Bins of Aj 

To study a potential effect of parton energy loss on the fragmentation properties of partons in 
more detail, we divide the data sample into classes of dijet imbalance. Four Aj bins are cho- 
sen, which split the data sample into approximately equal parts: < Aj < 0.13, 0.13 < Aj < 
0.24, 0.24 < Aj < 0.35, and 0.35 < Aj < 0.70. The fragmentation functions are reconstructed 
separately for leading and subleading jets in central events. Figure [5] shows the fragmentation 
functions in bins of increasing dijet imbalance {Aj) from left to right (top row), the correspond- 
ing ratio to the fluctuation convoluted pp reference (middle row) and the jet pj distribution for 
leading and subleading jet (bottom row). For both leading and subleading jets, independent 
of the dijet imbalance bin, the PbPb fragmentation functions closely resemble those of the pp 
reference. 



5 Summary 

The CMS detector has been used to study jet production in PbPb collisions at ^/s^= 2.76 TeV. 
No visible modification of the dijet azimuthal correlation due to the nuclear medium was ob- 
served while a strong increase in the fraction of highly unbalanced jets has been seen in central 
PbPb collisions compared with peripheral collisions and model calculations. This observation 
is consistent with a high degree of jet quenching in the produced matter. A large fraction of the 
momentum balance of the unbalanced jets is carried by low^-px particles at large radial distance 
to the jet axis. The hard component of fragmentation functions reconstructed in PbPb colli- 
sions for different event centrality and dijet imbalance Aj exhibit a universal behavior closely 
resembling jets of the same reconstructed energy fragmenting without having traversed a nu- 
clear medium, as seen in the comparison with pp collisions. These results provide qualitative 
constraints on the nature of the jet modification in PbPb collisions and quantitative input to 
models of the transport properties of the medium created in these collisions. 
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